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We examine Big Bang nucleosynthesis (BBN) in the case of inhomogenous neutrino degeneracy, in 
the limit where the fluctuations are sufficiently small on large length scales that the present-day ele- 
ment abundances are homogeneous. We consider two representive cases: degeneracy of the electron 
neutrino alone, and equal chemical potentials for all three neutrinos. We use a linear programming 
method to constrain an arbitrary distribution of the chemical potentials. For the current set of 

^\ , (highly-restrictive) limits on the primordial element abundances, homogeneous neutrino degeneracy 

^ ' ' barely changes the allowed range of the baryon-to-photon ratio rj. Inhomogeneous degeneracy al- 

lows for little change in the lower bound on rj, but the upper bound in this case can be as large as 

^^ ■ t; = 1.1 X 10"** (only i^e degeneracy) or 77 = 1.0 x 10"" (equal degeneracies for all three neutrinos). 

For the case of inhomogeneous neutrino degeneracy, we show that there is no BBN upper bound on 

CTS ' the neutrino energy density, which is bounded in this case only by limits from structure formation 

and the cosmic microwave background. 
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I. INTRODUCTION 

Although the standard model of Big Bang nucleosyn- 
thesis (BBN) is highly successful (for a recent discussion, 
see references KM) rnany variations on this model have 
been proposed [Bf. One of the most frequently investi- 
gated variations on the standard model is neutrino de- 
generacy, in which each type of neutrino is allowed to 
have a non-zero chemical potential |4| - [|ll| , and a num- 
ber of recent models have been proposed to produce a 
large lepton degeneracy ||l2| - Q. 

More recently, Dolgov and Pagel have suggested the 
possibility of inhomogeneous neutrino degeneracy [[l5| . 
Their model was proposed to explain the apparent dis- 
crepancy between various measurements of the primor- 
dial deuterium abundance in high-redshift Lyman-alpha 
clouds. Here we consider a more mundane possibility: 
that the neutrino chemical potential is inhomogeneous, 
but on much smaller scales. In particular, we assume 
that the amplitude of the inhomogeneities is small on 
length scales larger than the typical baryonic diffusion 
scales after nucleosynthesis, so that the element abun- 
dances are homogeneous today. We calculate the ele- 
ment abundances for this scenario and compare to obser- 
vational limits. Using a method similar to that in refer- 
ence pq| , we can simulate arbitrary distributions of the 
neutrino chemical potential, and so determine the upper 
and lower bounds on the baryon-to-photon ratio rj in this 
model. 

In the next section, we discuss our model for inhomoge- 
neous neutrino degeneracy and its physical consequences. 
In Section 3, we use a linear programming technique to 
calculate upper and lower bounds on rj in this model. Our 
conclusions are summarized in Section 4. We find that 
when the chemical potential is inhomogeneous, there are 
no BBN limits on the overall neutrino energy density; 
the only limits in this case come from other cosmologi- 
cal considerations such as structure formation |l^,|l^ or 
the CMB 18 1 . Not surprisingly, inhomogeneous neutrino 
degeneracy allows for a wider range of values for rj than 
does homogeneous degeneracy. 



II. INHOMOGENEOUS NEUTRINO 
DEGENERACY 

Consider first the case of homogeneous neutrino degen- 
eracy. In this case, each type of neutrino is characterized 
by a chemical potential /i^ {i — e, /i, t), which redshifts 
as the temperature, so it is convenient to define the con- 
stant quantity ^i = fii/Ti. In terms of ^i, the neutrino 
and antincutrino number densities are given by 
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while the total energy density of the neutrinos and an- 
tineutrinos is 



P= TT^T^ 
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Degeneracy of the electron neutrinos alters the n ^^ p 
weak rates relevant for BBN through the number densi- 
ties given in equations (||) and (g), while the change in 
the expansion rate due to the altered density in equation 
(0) affects BBN for degeneracy of any of the three types 
of neutrinos (see, for example, reference [[lO| for a more 
detailed discussion). 

What happens if this degeneracy is not homogeneous, 
as assumed in almost all previous work, but instead varies 
with position? Dolgov and Pagel |]15[ considered such a 
model in order to explain the discrepancy in observed 
deuterium abundances at high-redshift. In their model, 
^ varies on scales ~ 100 — 1000 Mpc, producing an ob- 
servable inhomogeneity in the present-day element abun- 
dances. We make the opposite assumption: we take the 
variation in ^ to be small on such large scales, and large 
on much smaller scales, so that elements are well-mixed 
before the present day, erasing any detectable inhomo- 
geneities. 

Although models have been proposed which produce 
inhomogeneities in ^ (see, for example, the discussion in 
reference [u9[), we will follow the example of reference 
|l5| and keep our discussion as general as possible. In 
general, one would expect a distribution of fluctuations 
in ^ over all length scales. However, since the neutrinos 
are relativistic, they will free-stream and erase any fluc- 
tuations on length scales smaller than the horizon at any 
given time. We will make only two assumptions concern- 
ing the fluctuations in ^: that the fluctuations are sig- 
nificant on large enough scales to avoid being erased by 
free-streaming, and that they become negligible on small 
enough scales that the resulting element distribution is 
homogenous today. 

The first of these conditions requires that the fluctu- 
ations are significant on scales larger than the horizon 
scale when the n <-> p reactions freeze out at T '^ 1 MeV. 
If this were not the case, then free-streaming would erase 
all of the fluctuations in ^ before BBN began. This hori- 
zon scale corresponds to a comoving length scale ^100 
pc today. 

The condition that the element distribution be homo- 
geneous today requires that the fluctuations in ^ decrease 
sufficiently quickly with length scale that they have no 
significant effect on nucleosynthesis on scales above the 
element diffusion length. Although no detailed studies of 
element diffusion have been performed in connection with 
inhomogeneous BBN scenarios, it seems safe to assume 
that complete mixing of the primordial elements will oc- 
cur on scales well within the nonlinear regime today, < 1 



Mpc. By requiring the fluctuations in ^ to be negligible 
above this scale, we can also ignore any constraints from 
CMB observations, which severely constrain models with 
fluctuations on larger scales |15|] . 

Given these assumptions, we can assume that BEN 
takes place in separate horizon volumes, with the value 
of ^ being homogeneous within each volume. At late 
times, the elements produced within each volume mix 
uniformly to produce the observed element abundances 
today. 

Note that with this set of assumptions, it is no longer 
meaningful to talk about the value of ^ for the neutrinos 
at the present. Since the neutrinos from different horizon 
volumes diffuse freely up to the scale of the horizon, the 
different thermal distributions with different values of ^ 
will combine to give a highly non-thermal neutrino dis- 
tribution at late times. Thus, in the inhomogeneous sce- 
nario, it is still possible to put constraints on the present 
value of Pi,, but it is meaningless to discuss limits on 
^, since the present neutrino distribution will be non- 
thermal and cannot be characterized by a single value of 
^. This effect is present at some level even at the time of 
nucleosynthesis. The neutrinos remain in thermal equi- 
librium down to a temperature T ~ 2 — 3 MeV, so that 
the neutrino distribution remains thermal down to this 
temperature, with a single unique value of ^ in each hori- 
zon volume. However, the n <-> p rates do not freeze out 
until T ^ 1 MeV, so that neutrino free-streaming after 
decoupling will tend to produce a somewhat non-thermal 
background as early as the beginning of nucleosynthesis. 
We have neglected this effect, which will be negligible in 
any case if the neutrino inhomogeneities are conflned to 
comoving scales > 100 pc. 

We also assume for simplicity that rj remains uniform 
in the presence of an inhomogeneous lepton distribution. 
This need not be the case if, for example, baryogenesis is 
related in some way to the lepton number 



III. THE EFFECT ON BIG BANG 
NUCLEOSYNTHESIS 

Given our discussion in the previous section, we as- 
sume that the distribution of each type of neutrino is 
homogeneous within a given horizon volume during nu- 
cleosynthesis, and characterized by a single degeneracy 
parameter ^i {i = e, fi, r). Different horizon volumes 
may have different values of ^i, so we characterize the 
distribution of ^i by a distribution function /(^i), which 
gives the probability that a given horizon volume has a 
value of ^i between ^i and £_i + d^i . 

What form should we choose for /(Ci)? In analogy 
with the distribution of primordial density perturbations 
(and in accordance with the central limit theorem) the 
most obvious choice is a Gaussian distribution. However, 
we can consider a more general case than this. Using lin- 
ear programming techniques like those in reference p6[. 



it is possible to analyze the general case of an arbitrary 
distribution /. Consider flrst the case where only Vg is 
degenerate, and suppose that /(fe) is an arbitrary dis- 
tribution. Then all of the element abundances will be 
functions of ^e (for fixed 77), and we can write, for a given 
nuclide A, 



Xa = 



XAi^e)f{^e)d^e 



(4) 



where A^(^e) is the mass fraction of A as a function of 
fe, and Xa is the mass fraction of A averaged over all 
space; after the matter is thoroughly mixed, Xa will be 
the final observed primordial element abundance. 

In order to test all possible distribution functions /, we 
can divide the range in ^e into discrete bins (not neces- 
sarily all of the same size) , and approximate the integral 
in equation (0) as a sum: 



A^ = ^A^,/,Ae 
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where the dependence of Xa and / on ^e is expressed 
through their dependence on the bin number j. 

For each of the elements of interest ^He, D, and ^Li) 
we have an upper and a lower observational bound. Thus, 
for each of these three elements, we can write down equa- 
tions of the form: 

Slower bound ^ / ^ ^jJj ^f^j ^ ^upper bound- V^j 



Furthermore, f{S,e) is normalized to unity, so 

If we now define 

then equations (o) and (^ become: 

^ lower bound ^ / ^ ^jPj ^ ^ upper bound; 
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If we put an upper and lower cutoff on these sums, so that 
we retain only a finite number of terms, then equations 
(y) and (|1Q| ) are in the form of the constraint equations 
in a linear programming problem, with the N indepen- 
dent variables being the Pj's. In reference |1^, the vari- 
able under consideration was 77 rather than ^, so that 
the final quantity which needed to be maximized or min- 
imized was the mean value of rj. In our case, we wish 
to determine, for a given value of 77, whether there is a 



solution to equations (^) and (|lO|). Since there are non- 
BBN limits on p^, we have chosen to take the quantity 
p'nl Pv as our objective function, where p'^ is the final 
mean total neutrino density in the degenerate case, and 
Pu is the neutrino density in the absence of degeneracy. 
(These densities include all three neutrinos and antineu- 
trinos). We then determine whether a solution exists to 
our constraint equations for a given value of 77, and scan 
through the allowed range of r\ until we reach an up- 
per and lower value of r\ for which a solution no longer 
exists. At these limiting values for 77, our linear program- 
ming routine gives the minimum possible value of p'^/ p^^ 
which we can compare to other constraints. 

We consider two representative cases of interest: first, 
the case where ^e 7^ and ^^^ = ^r = 0, which is equiva- 
lent to ^e ^ Ca" ^T' ^^^ ^^^ *^^^^ Ce = Cp = Cr ■ The latter 
is probably the most physically realistic case y,|l3|- Al- 
though we have discussed our linear programming proce- 
dure only for the case of v^ degeneracy, it generalizes in 
an obvious way for the case where £,e — S,^^ — £,t ■ We have 
not considered the most general possible case, in which all 
three degeneracy parameters vary independently. How- 
ever, as we shall see, arbitrary inhomogeneity in ^e alone 
allows absurdly large values of 77 to be compatible with 
BBN, so there is nothing further to be gained in consid- 
ering the most general case. 

We use for our limits on the element abundances the 
values in the recent review in reference |y|. For the pri- 
mordial helium-4 mass fraction, Yp, we take 



0.228 < Yp < 0.248. 



fll) 



The limits on the number ratios of deuterium and 
lithium-7 to hydrogen are: 



2.9 X 10"^ < D/H < 4.0 X 10"^ 



and 



1.3 X IQ-i" < ^Li/H < 2.0 X 10" 
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However, a BBN calculation with these limits alone yields 
no single value of 77 consistent with all three sets of limits. 
One can argue either that the theoretical uncertainties 
are large enough to account for this discrepancy [||, or 
that one of these sets of limits (most likely lithium) does 
not represent the true primordial abundance [y . We have 
chosen the former approach. Folding in the theoretical 
uncertainties in the BBN predictions from reference [gj, 
we take the following limits on D/H and ^Li/H: 



2 X 10"^ < D/H < 5 



X 10 



-5 



and 



1 X 10"^*^ < ^Li/H < 4 X 10" 



(14) 



(15) 



We have ignored the theoretical uncertainty in hclium-4 
because it represents a much smaller fractional change in 
Yp. 



We wish to emphasize that our general results are fairly 
insensitive to small changes in the limits quoted above. 
Since we are exploring a rather radical change to the 
standard model, we make no effort to perform an ultra- 
high-precision calculation. 

We used the procedure discussed above to determine 
the largest and smallest values of rj which are consistent 
with the limits on Yp, D/H, and ^Li/H in equations (|ll]), 
(p4|), and (15). Our mixing procedure requires the use of 
mass fractions, rather than ratios to hydrogen, so we have 
made this conversion in our calculation. 

Consider first the "standard model" with no degener- 
acy. For the limits quoted above, we obtain bounds on 
7? of 3.7 X 10^1" < 7y < 5.3 x 10^^". Now what hap- 
pens if we add a homogeneous neutrino degeneracy? We 
have calculated the bounds on 7/ for the case in which 
^e (only) can have an arbitrary value, and for the case 
where i,e — S,ti — £,t can be set to any desired value. 
For both cases, we find that the bounds on rj are al- 
most unchanged. (The lower and upper limits are both 
enlarged by less than 2%). While this might seem sur- 
prising in light of earlier similar calculations p(| , it is 
a consequence of the increasingly narrower limits on the 
primordial element abundances. With such sharp limits 
as those considered here, even a free variation in ^e or in 
S.e = S,fj, = S.T cannot significantly alter the limits on 77. 
(We could obtain a larger range in 77 by allowing Ve and 
either j/^ or Vr to vary independently, but we would still 
expect a narrower allowed range than in reference llG] 
because of the improved observational limits) . 

Now we proceed to the case of inhomogeneous degener- 
acy. As we have noted previously, there is no well-defined 
mean final value of ^i in this case, since the neutrinos 
mix at late times to produce a non-thermal distribution. 
However, the mean final value of p^, is still well-defined, 
so we can attempt to constrain it with BBN. Consider 
first the case of ^e ^ Cm'^t- In this case, all of the el- 
ement abundances go to zero in the limit of large ^e- 
Thus, if we take f{£,e) to have the form /(^g = 0) ~ 1, 
and /(^g = ^0) = /o <c 1, where ^0 is a sufficiently large 
value of ^e such that all of the element abundances are 
neglible, then as we take the limit where ^0 ~> 00, the el- 
ement abundances approach their values in the standard 
nondegenerate model, while pi,^ goes to infinity. Thus, 
in the case of inhomogeneous Ve degeneracy, there is no 
BBN limit on p^^. Of course, there are other cosmo- 
logical limits on p,y in this case, from the requirement 
that structure formation not be disrupted by the extra 
radiation |1Q,&^ and that the extra radiation not distort 
the CMB fluctuation spectrum |Q. Our argument also 
applies to the case where all three neutrinos have equal 
chemical potentials. 

There are still interesting limits to be placed on rj. To 
determine these limits, we calculated the BBN element 
abundances for a grid of values of £,. We took ^ in steps 
of A,^ =1.0 between ^ = —60 and ^ = 10. Wc embedded 
a smaller grid between ^ = —1.0 and £, — 1.0 in steps 
of A^ = 0.05. In calculating the element abundances for 



the degenerate case, we used the approximation given in 
reference |1C|] for the decrease in the neutrino tempera- 
ture at large ^. Although rough, this approximation is 
adequate for our purposes. 

For the case ^e 3> Cm'^t' ^^ ^^'^ acceptable solutions 
for rj in the range 



3.0 X 10"'" <?7 < 1.1 X 10"^ 
while for the case £,e — ^f^ ~ £,t, we have 

3.1 X 10"^" <?7 < 1.0 X IQ-^ 



(16) 



(17) 



The actual ^ values of the non-zero bins, along with the 
corresponding values for pj and p'^/pu are given in ta- 
bles I and II. Note that our linear programming method 
will always yield a final optimal distribution for the pj 's 
in which at most seven of the bins are non-zero (since 
equations (||) and (|^ correspond to a total of seven con- 
straint equations); effectively, this corresponds to a final 
distribution for /(^) which is a sum of at most seven 
delta functions (see references |0,Q for a more detailed 
discussion) . 

We see that allowing for a free distribution of the de- 
generacies significantly increases the upper bound on rj, 
particularly for the case of ^e ^ Cm'^t' t)ut decreases the 
lower bound only slightly. Furthermore, the minimum 
increase in the neutrino density needed to achieve these 
lower bounds is inconsistent with both structure forma- 
tion considerations [|0| and CMB observations ||l^. On 
the other hand, the value of p'^/p,^ needed to achieve the 
upper bounds on 77 is well within the regime allowed by 
both structure formation and the CMB. 



IV. DISCUSSION 

Our results indicate that, not surprisingly, the intro- 
duction of inhomogeneous neutrino degeneracy allows for 
a much wider range of rj within the constraints of BBN. 
Current limits on the primordial element abundances are 
so tight that even models with homogeneous degeneracy 
are tightly constrained. Similarly, inhomogeous neutrino 
degeneracy does not allow for a significant decrease in 77, 
and such models tend to give a neutrino energy density 
in conflict with other cosmological limits. It is quite im- 
pressive that even with the radical model discussed here, 
the limits on the primordial element abundances have be- 
come so tight that a significantly lower value of r/ cannot 
be achieved. On the other hand, inhomogeneous neutrino 
degeneracy can increase the upper bound on rj to quite 
large values: up to 77 = 1.0 x 10~^ for the case of equal 
degeneracies in all three neutrinos, and ij = 1.1 x 10~^ if 
only Ve is degenerate. These correspond to ilbh"^ = 0.036 
and 0.40, respectively. 

The distributions of ^ which produce these extreme 
values for 77 do not correspond to physically likely models. 
The usefulness of our linear programming calculation is 



that it allows us to establish upper and lower bounds 
on 77 for arbitrary distributions of ^, while at the same 
time giving the smallest value for p'^/ Pu corresponding 
to a given value of 77. Any other distribution /(^) is 
guaranteed to give values of 77 which lie inside of our 
bounds. 

This work does not exhaust the possible models with 
spatially-varying ^. It is possible to use our methodol- 
ogy to investigate models in which two or three neutrino 
degeneracies are independent. In addition, if baryogene- 
sis is related to the lepton number in some way |0,n3, 
then one would expect a correlation between 77 and ^ at 
each point in space. Given a specification for this cor- 
relation, such models could also be examined within the 
framework we have outlined here. Less general but more 
physically realistic distributions for /(^) (e.g. a Gaussian 
distribution) could also be considered. 
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TABLE II. The nonzero values of ^_,, along with the 
corresponding pj and total p'^lpv for the case where 

e = Ce = ^M = e.- 



■n 


i, 


Pi 


p'uIpv 


3.1 X 10""' 


-0.7 

9 

10 


0.52 
0.30 
0.18 


103 


1.0 X 10"^ 


-25 

-1 

2 

3 


9.3 X 10"^ 
0.37 
0.45 
0.18 


3.2 



TABLE I. The nonzero values of ^j, along with the corre- 
sponding Pj and total p'^/ Pv for the case where ^ = .^e, and 

^M = Cr = 0. 



n 


6 


Pi 


Pulpy 


3.0 X lO"^'^ 


-0.8 
8 
9 


0.49 
0.03 
0.48 


31 


1.1 X 10"* 


-37 

-2 

3 

4 


7.3 X 10-^ 
0.23 
0.30 
0.47 


4.3 



